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1. Introduction. 

The inclusive cross sections of the leading neutrons in deep inelastic scattering with protons (DIS) attracts 
interest since this process allows us to extract the F2 deep inelastic structure function for pions {F2) (see 
Fig. and Refs. [1] ). This process has one rapidity gap since no particles are produced in the region of 
rapidity between the neutron and the bunch of secondary particles with mass {Mx in Fig. |l|). For long 
time it has been known that the cross section of such processes has to be multiplied by factor S'^ which is 
called survival probability [2-4] . This factor stems from the possible interaction of the constituents of the 
projectile with the target that should be forbidden to preserve the gap. In other words, the constituent of 
projectile could interact with the target in initial or final state suppressing the cross section of such process 
(see Fig. ^ that illustrates such interactions). 

The two decades experience in calculation of the survival probability ( see Refs. [5-11]) has shown 
that the predictions for survival probability could differ in the orders of magnitude (from S'^ = 0.9 to 0.01) 
depending on the model for strong interaction at high energy. 

In this paper we develop the theoretical approach for survival probability in DIS based on high density 
QCD [12-17]. For DIS at high energy QCD predicts that the dynamics of quarks and gluons (partons) can 
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Figure 1: The generic picture for leading neutron production in DIS: the Born approximation.) 

be described in terms of parton saturation* [12-14]. The new state of the matter: Color Gluon Condensate 

(CGC) [14,15] will be produced in which the amplitudes are dominated not by quantum fluctuations, but by 
the configurations of classical field containing large, ~ 1/05 numbers of gluons. This state is characterized 
by large density but small QCD coupling as- This feature results in the solid theoretical approach based 
on non-linear equations [12,13,15-17]. It worthwhile mentioning that this theoretical approach has reached 
the most developed stage for DIS with which we are dealing in this paper. 




The main result of the paper is the equation for 
the inclusive production of leading neutron. The basic 

theoretical ideas how to approach the processes that 
have the rapidity gap, have been formulated in Ref. 
[18]. We explore these ideas to obtain the equation. 
We found the solution to the new equation and show 
that the interaction in initial and final sate will lead 



to the cross section that falls down at high energy. In 

Figure 2: Leading neutron production in DIS: general ^^^^^ ^^^^^^ survival probability turns out to be 

source of shadowing corrections , , . , 

small at high energy. 

The paper is organized as follows. In the next section we will derive the equation and discuss the high 
energy asymptotic behaviour of the solution. Section 3 is devoted to the initial condition to the equation, 
while in section 4 we obtain the numerical solution to the equation. In conclusion we summarize the main 
results. 



2. The equation 

2.1 Derivation of the equation in the dipole approach 

In inclusive leading neutron production (ILNP) the produced partons could interact between themselves 

as well as with the target in the final state, unlike in the case of total cross section for which we have 

the non-linear Balitsky-Kovchegov equation [16, 17] . In principle, the interaction in the final state means 

that the cross section of the process does not depend only on the structure of the wave function of the 

colliding particle as it is for the total cross section. However, in Ref. [18] (see also Ref. [19]) it is shown 

on the example of diffraction production, that in the processes which are inclusive, the different type of 

interactions in the final state cancel each other. In this section we argue that ILNP is also belongs to the 

class of such inclusive processes as diffraction production. 

* Parton saturation means that at high energy the density of partons (actually mostly gluons) reaches the maximum value. 
It is instructive to notice that in CGC the QCD coupling is small (as 'ti 1). 
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In ILNP we can see three typical moments of time^ : at x~ = — oo we have fast virtual photon (or, 
better to say, a system of partons in the coherent state of the virtual photon) and the target, at x~ = 
these partons interact with the target and the produce partons that propagate to x~ = +00 where the 
detectors are placed which measure this system of partons. Actually in our case we do not measure these 
partons which means that we sum over all possible final states with only one restriction that their total 
mass is equal to Mx(see Fig. ^). The same time structure we have in complex conjugated amplitude in 
which we denote times as x~ (see Fig. ^). 




Figure 3: The time structure of the inclusive production of leading neutrons. The low panel shows our notation for 
the process shown in the upper panel. The straight solid lines denote the initial dipole in the low panel. 



Fig. ^ shows the particular contribution to the cross section of interest. From x^ = —00 to x~ = the 
wave function of the fast dipole consists of system of two quark-antiquark pairs and two gluons, at x~ = 
this system interacts with the pion which destroys its coherence. As a result the components of the system 
(two gluons and two quark-antiquark pairs ) are produced. The components interact in the time interval 
< x~ < +00. The gluons can be emitted and absorbed in the final state. Because of this we cannot 
use optical theorem that reduces the cross section to the imaginary part of the elastic amplitude. The key 
observation is that absorption and emission in the final state ( for times > and x^ > 0) cancel each 
other. This cancellation was first proven in Ref. [21]. It is found that we have two types of cancellations 
shown in Fig. ^. The first type is cancellation of the gluon emitted in the final sate and which can be 
caught by the detector (see the first diagram in Fig. ^A) with gluons emitted and absorbed in the final 
state in the amplitude and the conjugated amplitude (see the second and third diagrams in Fig. Q-A). The 
second type of cancellation is shown in Fig. ^B. The produced gluon from the initial wave function, that 
exists from x~ = —00 to x~ = +00, cancels by the gluon from the wave function that has been absorbed 
in the final state and, therefore cannot be measured by the detector. Notice that the gluon in the first 
diagram of Fig. Q-B has been emitted in the final state in the conjugated amplitude. In Fig. ^ we denote 
the gluon by quark and antiquark lines since we consider our process at large Nc where Nc is the number 
of colours. 

Armed with these cancellations we can show that in our process the emission and absorption in the 



^We will consider all processes in the light-cone frame and will use the light-cone perturbative theory (see Ref. [20]) 
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Figure 4: Two types of cancellations in the final state 



final state does not contribute to the cross section. To illustrate this point we consider the example of 
Fig. ^. One can see from Fig. ^ that emission of the upper gluon in Fig. |5|-1 is canceled by diagram of 
Fig. ^2 which describes the emission and absorption of the upper gluon in the final state (actually for 
cancellation we need to add the diagram of Fig. ^2 type for the conjugated amplitude). The absorption 
of the low gluon in Fig. |5|-1 is canceled by the diagram of Fig. |-3. This dia gram describes the emission 
of the gluon from the initial state in the amplitude but the emission of gluon from the final state in the 
complex conjugated amplitude. The cancellation of emission and absorption in the final state means that 
we can used the optical theorem for interaction of the dipole with virtual pion (see Fig. |l]). It is worthwhile 
mentioning that this cancellation works in any complex diagrams which differs only by the gluon shown in 
Fig. I [21]. 

The formal proof of the discussed cancellation can be done using the method of mathematical induction. 
Let us assume that for emission of n-gluons we have the cancellation and they can be characterized by the 
initial wave function of colorless dipole. The emission of n + 1 gluon by one dipole can be described by 
the sum of the diagrams of Fig. One can see that diagrams of the set A in Fig. ^ cancels as well as 
the diagrams of set B , due to the cancellation give by Fig. ^ The only diagram that remains and gives 
the contribution is the diagram of Fig. |6|-C. Since the initial condition for N"^ contains only one dipole we 
conclude that we can describe this processes by the initial wave function of dipoles. 

Now we can derive the equation based on the same ideas as have been used in the derivation of 
Balitsky-Kovchegov (BK) equation, since the ILNP stems from the structure on the initial wave functions 
of colliding particles. Let us denote by N"^ (x, b, Y, Yg) the imaginary part of the amplitude shown in Fig. |l|. 
Yg denotes the rapidity gap between the neutron and the slowest hadron among the produced particles 
(see Fig. I). 

A long before the interaction the incoming dipole decays into two dipole each of them could scatter 
separately and produced the leading neutron. These two dipoles can interact simultaneously producing 
leading neutrons. However, we have another process which does not produced any particles that fill the 
rapidity gap. One dipole scatters elastically while the second produces the leading neutron. Fig. |^ shows 
these three ways of interaction. All these terms we can see in Fig. |^ The first term shows the separate 
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Figure 5: The example of cancellation of emission and absorption in the final state for the contribution of the 
diagram of Fig. ^ to the cross section. 



interaction of each dipole while the second term describe the simultaneous interaction of two dipoles. The 
third and the forth terms describe the ILNP by one dipole while the second dipole undergoes the elastic 
scattering. The analytical form of this equation looks as follows: 
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Figure 6: The emission of one gluon by the dipole. The diagrams with the gluon emitted and absorbed by the same 
quark (antiquark) are not shown but they have the same pattern. 
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Figure 7: The graphic form of the equation. The vertical dotted hnes denote the time of interaction in the ampUtude 
or in complex conjugated amplitude. If such line crosses the blob it means that the interaction occurs. The wave 
lines denote the infinite time and if such line goes through the blob it means that the blob describe the production 
of particles. 



dY ^ ' ' 



^d^^2^^ I^N^ (^xo2,b - lxi2;Y,Y/\ + ( xi2,b - lxo2;Y,Yg 



27r J Xq2 

(xu, b; Y, Yg] + { xo2, b - ^xis; Y, Y„ ] { xi2, b - ^xo2; Y, Y. 



- 2N^ (^xo2,b-^xi2-,Y,Yg^ N (^xi2,b-^xo2;Y'^ - 2N (^xo2,b - ^Xi2-,Y^ (^xi2,b - ^xo2;Y,Yg^^ 
where N (xio, b, Y) is the solution of Balitsky-Kovchegov equation. 



2.2 BFKL Pomeron calculus and inclusive production of leading neutron 

As it is well known the non-linear BK equation for the elastic amplitude corresponds to summation of 'fan' 
diagram in the framework of the BFKL Pomeron calculus [12,22]. The equation for the ILNP stems from 
the AGK cutting rules [23]. These rules has been proven in QCD for all processes in which we do not have 
emission of the gluon from the triple BFKL Pomeron vertex (see Refs. [24-29]). Our process belongs to 
this class and we can safely apply the AGK cutting rules. The contribution to the imaginary part of the 
elastic scattering amplitude can be viewed as sum to three different processes which are shown in Fig. ^. 
If the rapidity Y' is the position of the triple Pomeron vertex the simplest 'fan' diagram generates three 
different processes: in the first one (double cut in Fig. ^ the state with double density of particle (gluons) 
is produced while in the second process (single cut in Fig. |8|) the density of the particle is the same as in 
the single Pomeron exchange; the third process is the diffraction production in which we do not produce 
a particle in this rapidity gap. The AGK cutting rules establish the relations between these processes and 
they are shown in Fig. ^-A. 



Cut P(.inier()n 



Diffract I ve cut 



Fig. g-a Fig. g-b 

Figure 8: The processes that corresponds to triple BFKL Pomeron contribution to the total cross section. 
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The difference between coefficients that 
are shown in Fig. |^A and the standard ones 
is related to the fact that we use for the con- 
tribution of the cut Pomeron the contribu- 
tion to the total cross section (see Fig. §-a). 
instead of the contribution to the imaginary 
part of the amplitude. As we know from the 
unitarity constraint: 2Imj4 = atot- 

Having Fig. ^-A in mind we can see that 
the leading neutrons can be produced from 
the first and second contributions while the 
third one corresponds to the diffraction pro- 
duction and it contributes only to the spec- 
trum of leading protons (see Fig. ^B). Re- 
placing Pomerons by the the sum of 'fan' 
diagrams we obtain the equation given by Eq. (^). The direct relation between dipole language in QCD 
and the BFKL Pomeron calculus for the inclusive processes have been noticed in Ref. [18]. The powerful 
theorem on cancellation of the interaction in the final state in the Pomeroin language are hidden in the 
assumption that only Pomerons and their interaction contribute to the inclusive processes. Therefore, the 
dipole consideration can be considered as the theoretical argument for the Pomeron calculus. 



Figure 9: The relation between the AGK cutting rules and con- 
tribution to the inclusive production of leading neutrons. The 
helix lines denote gluons, pious are shown by the dashed lines. It 
should be noticed that in Fig. between two exchanges of pions 
the thick line denotes the neutron. 



2.3 Solution at ultra high energy 



We need to find an initial condition to Eq. {2A) before searching for the solution. However, the experience 
with BK equation shows that the asymptotical behaviour at high energies does not depend on the initial 
solution [30]. Using the approach of Ref. [30] we can find the solution noticing that the solution to BK 
equation approaches unity at ultra high energies. Substituting = 1 in Eq. ( |2.1| ) we obtain the following 
equation for A^'^ 
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^-02 -^12 



- [xi2,b-Y,Yg) + (xo2,b-]^xi2-,Y,Y^ (xi2,b-]^XQ2-,Y,Yg 



Assuming that N'^ falls down at high energy we can neglect the (Ntt)'^ - term and one can see that Eq. (|2^ 
can be re-written in the form: 



^A^" {xio,b;Y,Yg) 
dY 



(2.3) 



1 ^ \ , _ / r 1 



N 



^ (xi2, b; Y, Yg)) -2 {xoi;Y, Yg) } 



This equation can be solved using the double Mellin transform for N, namely, 

f-e+ioo J fe+ioo jf 

- L ^.L (2.4) 

where ^ = In {x'qi ) . It should be stressed that we need to consider x\q (Y, 6) S> 1 since only in this 
region we can replace by unity (Qg is the saturation momentum for A^'^). 

Indeed, for we obtain the following equation 

Lon^{f,b;LO,Yg) = -as{x{f) n^f,b;co,Yg) +4 ' (2.5) 

where x if) = 2^(1) — ip{f) — — f) is the BFKL kernel [38] in /-representation is di- gamma 
function , see formulae 8.360 - 8.369 in Ref. [39]). 

The solution to Eq. ( |2.5| ) is 

fi-{f,b;u,Yg) = C{Yg,b)e^^s+-^^^''^'^^f'^ (2.6) 



Taking Mellin integral of Eq. ( |2.4| ) by the steepest decent method we obtain the following equations for 
the saddle points: 

if +ix «^0; (2.7, 

Taking into account that x if) '^'^^> —2 In / we obtain the solution: 

N"" {tb;Y,Yg) = C7(yg,6) e-^'^s^^^-''^^) (2.8) 
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This solution demonstrates that the influence of the interaction in initial and final states are so essential 
that they lead to decrease of the cross section of the inclusive production of the leading neutrons. This 
makes the extraction of pion deep inelastic structure from this experiment problematic if not impossible. 
However, for more practical conclusions it is necessary to solve Eq. ( |2.lD in more realistic kinematic region. 
For this goal we need to know the initial condition that we are going to discuss in the next section. 



3. Initial conditions 



The problem of the initial conditions actually includes two different issues. The first is to find the expression 
for the Born Approximation (see Fig. |l|) for the leading neutron production in DIS for the initial energy. 
The second one relates to the explicit form of shadowing correction at the same low energy. Restricting 
ourselves by calculation the survival probability we do not need to know the details of Born Approximation 
which can be found in Ref. [1], since they cancels in the expression for the survival probability, namely, 

^2 c?(Tj„cZ (exact with shadowing corrections) / dcrj„c« (Born Approximation of Fig. |l]) (3 1) 

dYnCpPn,T I dYnd'^Pn,T 

where Yn and pn,T are rapidity and transverse momentum of produced neutron. 

The Born approximation at low energies looks as it is shown in Fig. [l^. The expression for this diagram 
takes the following form (see Ref. [1] ): 



z 



dzd^qr 167r2 (^rnl + q^) Q 



where aem is the fine structure constant. In the rest frame of proton the longitudinal momentum, 
carried by pion , is equal to ql = {I — z)m,^ j ^ and 



Xn = Y~c ^ Q2 _^ j^^2 ' ^9 = -ln(l - z)\ = a {m\ - q^) pion Regge trajectory (3.3) 

For initial condition we need to take F'^ (xj^ — )• 1, Q^) at — ?• 1. For calculation of survival probability 
(see Eq. ( |3.lD ) we do not need to know the value of F'^ and even its x-,^ dependence. The only ingredient 
that we need, has not appeared in Eq. (p.2|). The total cross section corresponds the 7*7r amplitude that 
is taken at momentum transferred from virtual photon to virtual photon (A in Fig. [I0| ) to be equal to 
zero. For calculation of the survival probability we need to know the A dependence of this amplitude. We 
explain this fact a little bit below. The expression for such an amplitude takes the following form 

N^{A;Q;Y = Yg,Yg) = (3.4) 



7QT ■ 
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We do not know the dependence of F'^ on A. In the additive quark model (see Fig. ^) it is natural to 
assume that F"" {x^;Q^;A) = F"" (x^;^^) x (A^) where is the electro -magnetic form factor of 
pion. On the other hand the main contribution for scattering at low energy gives the /^-resonance which 
contribution has no A-dependence. 

Integrating Eq. ( |3.4[ ) over qx and going to impact 
parameter representation 



7T /q- zr-cr\ 



N-{b;Y = Yg,Yg) = l^e'^^-'N-{A) (3.5) 
we obtain that 

Figure 10: Born Approximation for scattering ^ f ^f^_ ^2 /^a (3 g) 

art low energies. The dashed line describes pion, the J ^ J ^ ' 

solid arrowed line corresponds to quark(antiquark) 

where factor (Yg) includes the dependence of Ygi^z) and does not contribute to the calculation of the 
survival probability (see Eq. (|3.1|)). In Eq. (|3.6|) G^^ (b) is the impact parameter image of Gt^ (A) while 



^^^) = I / |^e--Vl-.)<^^^^%^ = J- / 4rf.TJi(.T5) %^(l-.)<^^ 
do J [2TTy + qf, 2?: J + qj. 

(3.7) 

We use (A) in the standard form 

°-<^' = rrhM 

where /i^ = 0.6GeV^ is taken from the measured value of the electro-magnetic radius of pion:i?^ = 
0.66ib0.01/m [31]. For Gp^+„ form factor there exists a variety of different experimental (phenomenological 
) information (see Ref. [1] for details). We take this form factor to be proportional to electro- magnetic 
form factor of proton in the spirit of the additive quark model that we pictured in Fig. IC. It takes a form 

Gpn+n (Qt) = 9p-K+n Gp {Qt) = " 2 ('^•^) 

(1 + Qt/lJ'l) 

with fil = 0.72Gey2 which corresponds to Rp = 0.862 ± 0.012/m. In Ref. [1] is suggested to replace 

Gp.+„(gT)e<^^^)^ .^1,2 ^ith J- = ^ + a;y, (3.10) 

This expression describes correctly the small qx behaviour which is the most essential feature of the pion 
exchange. 

Using Eq. ( 0.10 ) we obtain 

„2 

T (b) = 2 '5 2 ("^-^1 ("^-^) - W ^1 if'effb) ) (3.11) 



- 10 - 



We will use Eq. (|3.10| ) in our estimates. 

Now we need to specify the initial condition for BK equation. As it will be seen below we cannot 
calculate the survival probability without introducing the impact parameter dependence of the scattering 
amplitude. It is known, that BK equation has a problem with taking into account the large 6-dependence of 
the scattering amplitude [37] and should be modified at large values of b. We suggest a different approach 
in the spirit of the additive quark model, assuming that we have two different scales inside the proton: 
the size of the proton and the size of the constituent quark which is much smaller that the size of proton. 
Since the saturation scaler is larger that iZ-^constituent quark integrate all amplitude in the BK- 

equations over impact parameter. In this picture the entire 6-dependence will be concentrated in initial 
conditions and will be determined by the form factor of proton. In this model we view a proton in the 
same way as tritium in the Glauber approach (see Ref. [32]). Such an approach could be relevant only if we 
know that the radius of interaction between dipole and constituent quark does not increase with energy. 
In high energy phenomenology based on soft Pomeron approach, the increase of the radius of interaction 
R with energy looks as follows 

i?2 =i?^ + apln(s) < rI (3.12) 

where Rp is the size of the proton and Rq is the size of the constituent quark, a'p is the slope of the 
Pomeron trajectory. Fortunately, some models have recently been discussed with very small value of 
a'jp ^ O.OlGeV'"^ [33] and Eq. ( |3.12| ) could be valid at high energies. 



We will fix the initial condition using results from Ref. [34] in which a perfect fit to HERA DIS data 
was made in framework of BK equation. We use the same initial conditions as in this paper, namely,"'' 

GBWmodeh iV (xio, F = F,;) = ctq (l - exp (- (xgi Qg,) ^4) ) ; (3.13) 
MV formula: N (xio,r = Yg) = ao (^1 - exp (x^^ Ql,y In ( ^^^^^^^ + /4^) ; (3.14) 

but we introduce the b dependence considering that initial saturation momentum Qq^ depends on b in 
the following way: 

Ql Qlsib) = QlsGpib) = QlmbKiimb) (3.15) 

where Gp (5) is the b image of Eq. (3^). Therefore, Eq. ( p. 13 ) and Eq. ( |3.14| ) transform into the following 
expressions 



N{xw,Y = Yg) = J <fbN{xio,Y = Yg-b) = < 



{ JdH{l-exp{-{xlQlib)r/i)^ 

(3.16) 

f£b{l - exp (- (xg, Ql {b)y In + e) /a) ) 



* Golec-Biernat and Wusthoff (GBW) model being very simple, describe DIS data [35] while McLerran-Venugopalan formula 
is the correct initial condition in Colour Gluon Condensate approach [14] 
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Considering Eq. (|3.13D and Eq. ( |3.16|) at small xqi (xqi Qq^ ^ 1) we see that 



2 ^7 



d25(Q2^ (6))' 



(3.17) 



This equation gives us the value of Qos as function of fio and m. In our solution we took two values 
for Ql^ in Eq. first one, considering m being the same as in the electromagnetic form factor of 



proton = 0.72GeV^; and the second one, taking Qq^ in Eq. ( 3.13| ) and Eq. (|3.14|) being equal to Qq^ 



in Eq. ( 3.13| )(Eq. ( |3.14D ) but changing m to satisfy Eq. ( 3.17 ). Since it turns out that 7 is close to 1, we 
take as a first try Eq. ( |3.15 ) with 7 = 1. 

In this case we have simple initial condition for BK equation in momentum representation for Eq. ( |3.13 ) 



Nik,Y = Yg,b) 



xoiiixoi Jo (bxoi) 



N{xoi,Y = Yg,b) 



''01 



^ro {k'/Ql (b)) (3.18) 



where To is incomplete Euler gamma function (see formulae 8.25 in Ref. [39]). In the case of the initial 
condition of Eq. ( |3.14 ) we did not find a simple analytical form and perform the integral of Eq. ( |3.18D 
numerically. 

The easiest way to get the initial condition for A^'^ is to write them in coordinate representation. 
Assuming that at large -F^ = (2/3)^2 due to the quark counting rules [36] we can obtain the initial 
condition for N'^ in the form of Eq. ( |3.13D and Eq. ( |3.14| ) in which we replace Qos(^) (2/3)Qos(^)- 
other words the initial conditions for A^'^ have the form 



(3.19) 



f (l-exp(-|(x^,gL(6))V4) 



N^^ {xio,Y = Yg;b) 

Ng{Yg) G.p (b) I (Pb' (b-b'^ r2 (6') < 

^ (1 - exp (-1 {xl, Ql ib)yin + ej /4 

However, we need to multiply these initial values of N'^^ by the survival probability since the interaction 
in the initial state can fill the rapidity gap for the leading neutron production. 



In the case of Eq. ( 3.15| ) the survival probability is equal to 

For GBW model: = exp (^-^a^oi Qos Gp {b)^ ; 

For MV formula: = exp ( -^xgi Ql, Gp (b) In ( ^ 

V 2 V^io^i 



+ e 



QCD 



(3.20) 



since due to unitarity constraint such S'^ corresponds to probability not to have any inelastic interaction 
which can spoil the rapidity gap [2,4]. Therefore the initial condition for N'^ in the coordinate representation 
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has a form 

?2 AJ-TT 



For GBW model: = A^^ {Yg) (l - exp (^-^ (x^, Qg^ (6))^ /4^ 

X exp (^-ix^i Ql Gp (6)) ^ d^ft' (b - b') {b') (3.21) 

For MV formula: = Ng (Yg) ( 1 - exp ( -| (xgi Q^, (fe))^ln ( ^ + e ) /4 ) ) 

V V \xioAqcd J J ) 

xexp(^-ix^iQ2,Gp(6)lii(^^-i-— + e))y" (6') (3.22) 

Going to momentum representation we have the initial condition in the form: 

N"" {k,b,Y = Yg,Yg) = 2iT I xoidxoiJo{xoik) ^ ^ — ^ 



^01 



Factor A'^ (1^) is not important for the calculation of the survival probability if we neg lect (iV'')^-term 
in Eq. ( p.l| ) (see below). However, we need it for the estimates of the accuracy with which we can neglect 
this term. It is equal to [1] 



IGvr 

4. Numerical solution 



Ar^(y^) = ^^e-^« « 2.2e-^« (3.23) 



In this section we will discuss numerical solutions of two equation: BK equation and Eq. (2.1), in momentum 
representation 



N{xm.b,Y) = xl^ J kdkJoikxQi) N {k,b,Y); 

{xo,,b, Y, Yg) = xl, [ kdkJo (kxoi) {k, b, Y, Yg) ; (4.1) 



In this representation BK equation looks as follows 
dN{k,b,Y) 

Wy 

d^k\ 

as < 



(4.2) 



/ ,2 (N{k',b,Y) ^ N{k,b,Y)) -NHkAY)}; 

[k-k'j^ k'l + [k-k'j^ J 

This equation we solve with the initial condition given by Eq. ( 3.161 ). 



- 13 - 



After finding the solution to Eq. ( |4.2| ) we solve the following equation: 

{k,b,Y,Yg) _ 



dY 



(4.3) 



as < 



^ ^ N 2 

k-k'^ 



N"" (k',b,Y,Yg 



i^) ^; Yg] 



'iN{k,b,Y) N'^ {k,b,Y,Yg 



k'l + (k-k' 

with the initial condition given by Eq. ( [3.23 ) 

Eq. ( [4.3| ) differs from Eq. (2.1) in the momentum representation by the (iV^)^-term which we neglected 
in this equation. Indeed, A^'^ in Born approximation is much smaller than N both experimentally and 
theoretically, since it describes a specific configuration that contribute to the inclusive cross section given 
by N. This configuration is suppressed by factor e~^^ . The second argument stems from the solution at 
high energies (see section 2.2) which shows that {N'^)^ turns out to be much smaller than AN N'^. We 
proceed with solution of Eq. (4.3) but after finding the solution to this equation we will check that the 
(A^'^)^-term is small. 

Using the solution to Eq. ( |4.3| ) we can find S"^ for dipole scattering which is equal to 

= J(fbN^{k,b,Y,Yg) / f£bN^^{k,b,Y,Yg) (4.4) 

where A^'^ is the solution to Eq. ( [4.3| ) while N'^^ is given by Eq. ( [3.21 ). 

One can see that for self-consistent calculations we need to estimate the contribution of the Born 
Approximation using Eq. ( |3.6D as the initial condition, solving the linear BFKL equation. 

{k,b,Y,Yg) _ 



dY 



(4.5) 



as < 



(fk\ 



k-k' 



AT- {k',b,Y,Yg) 



k' 



kf + (k-k 



■N^ {k,b,Y,Yg 




Therefore, we need to solve numerically three equations: Eq. ( |4.2| ), Eq. ( ^.3] ) and Eq. (|47 
that using a new variable for GBW initial condition (see Eq. ( p.l3D ) 

k(6) = k/Qos{b) 

the initial conditions for these three equations can be written in the following form 



Eq. (4.2) 



Eq. (4.3) 



Eq. 



N{k,Y = Yg,b) = -To{K^{b)) ; 

Nik,Y = Yg,b) = Ng (Yg) [ <f b' (b " b') {b') 



Af" {k,Y = Yg,b) 



r \ 

NgjYg) 



e 6 



1 2 

e~2'' Jo(Kr) 



We notice 
(4.6) 

(4.7) 
(4.8) 



d'b'Gjb-b'] T'[b') -T,[-K'{b)]- (4.9) 
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Since Eq. (|4.3|) and Eq. (^^) are linear equations with respect to N'" we can find the solution to them 
with the simplified initial conditions 



Eq. (H 
Eq. (P 



{k,Y = Yg,b) 
ik,Y = Yg,b) 



r \ 
2 

' 



e 6 



e 2^ 



Jo(Kr) 



(6) 



(4.10) 
(4.11) 



and the solutions with the initial conditions of Eq. (4.8) and Eq. ( |4.9| ) can be obtained as 



N^{k,Y;b) 



Ql (Xg,b) 



J dH' (b - 6') {b') N"" {k (6) , Y) (4.12) 



where N (k (6) , Y) is the solution to Eq. (O) or to Eq. ( [4.5| ) with the initial conditions given by Eq. ( [I.IOD 
or Eq. (ITTll ). 

All equations that we are discussing here are conformal invariant and can be re-written in variables 
K and k' in stead of k and k' . Hence actually we need to solve the system of three equations (Eq. ( |4.2| ), 
Eq. (|1|) and Eq. (U)) with the initial conditions of Eq. {^), Eq. (|4.10D and Eq. (|Tl|) to find N {k,Y) 
and A^'^ {k,Y). In spite of simplicity of Eq. ( [4.7D .Eq. ( [4.10D and Eq. ( 4.11 ) they have clear shortcoming 
since cannot reproduce a correct behaviour accordingly to operator product expansion at large values of 
k'^. Indeed, they fall down exponentially instead of 1/k'^. Therefore, we can trust these equations only at 
k'^ < Q1- Trying to find a compromise between simplicity and rigorousness we chose the initial condition 
for N'^ in the form 



3 \xqiKqcd , 

which is the expansion of Eq. ( 3.19| ) for McLerran-Venugopalan formula at small values of xqi- We can 
only treat xqi < 1/Qs with this initial condition but since Qs for low energy in DIS with a pion is small we 
believe that we can use Eq. ( f4.13| ) as the first approximation. In our equation for N'" the solution of BK 
equation is essential only for xqi < Qs and, therefore, we can use the GBW model for the initial condition. 
The same situation in the initial condition for N'^^. Finally, we use 



In 



(4.13) 



BFKL equation: 
our equation: 



2^01 Qos 111 
-Xqi Qqs In 



1 



1 

Xoi ^QCD 



exp 



^01 Qos/^ 



These equations take the following form in the momentum representation: 



BK equation: 
BFKL equation: 
our equation: 



1 



-Fo (k'/Qi) ; 
2 Qos/^'^'i 



1 



1/a 



k^ 

hs 



exp 



k^/Q 



Os 



(4.14) 
(4.15) 

(4.16) 
(4.17) 
(4.18) 
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The expression in Eq. ( [4.18| ) follows from the following simple calculations: 

1 



1 t-^ I 1 \ Fourrier transform 1 
ln(l/a;oi) > 



Fourrier transform / I \ \ 



and from formula 8.357 of Ref. [39]. 

We need to find the pion deep inelastic structure function, using N'^ and A'^g^, for calculating the value 
of the survival probability. In the dipole approach the cross section for the virtual photon with the pion 
can be expressed through the amplitude of the dipole-pion interaction in the following way for massless 
quarks: 

cj^.^(Q,y) = y"rf^xoi|^'(Q;xoi)pcT^iipole.7r(a;oi,^) = J (fxoi\^ (Q; xoi)\'^ J (fbN^ {xqi, Y;b) 

(4.20) 

where [40] 

1^ (Q; xoi) \' = Y.^} I ^'i'' + (1 - ^)'] Q'^i (^^oi) (4-21) 

^ / 

where -^z is the fraction of the electrical charge aem carried by quark of flavour f;Qis the photon virtuality 
and = Q'^z{l — z) ; z is the fraction of photon energy carried by the quark an Nc is the number of 
colours. Using Eq. Eq. (|^) can be re-written in the form 

'^rniQ,Y) = / kdk^^^,{Q,k) I dHN^{k,Y;b) (4.22) 

where 



$2 J r = ^ 
^ \ k 



) = j ^'^01 / - ^) + (1 - ^)'] ^01 (kxoi) Kl {Q xqi) (4.23) 

1 - k2)V1 + 4^2 + 8^2(1 + K^) ArcCsch (2k) 



/ dzz{l- z)[z' + zY]- 



(1 +4k2)V1 + 4k2 



where k = Q/k = T^J z{l — z). 

Using N"" {k, Y) , Eq. (|]2|) and Eq. (g!23|) the final expression for the survival probability (see 
Eq. ( [4.4D ) has the following form 



Jd-^b JdWG^(b-b')T^{b') I 



^QCD/Qosib) 



KdK {Q/iQosib) iV- (Eq. (|;5D ; ^,Y) 



In Eq. ( 4.24 ) we use the infrared cuttoff k = Aqcd since for smaller k we cannot use Eq. ( p^ for 
the dipole amplitude. 
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ffsY= 0.4 




»sY=2 




»sY= 5.2 






ffjY=0 




ffsY= 0.4 




itsY=2 




ffsY= 5.2 




OsY=0 
QjY= 0.4 
OsY=2 
asY= 5.2 
QsY= 10 
QjY= 20 



Fig. 11-b 



Solutions to equations at different rapidities: 



Fig. |ll|-a shows tlie solution to the Balitsky- 
Kovchegov equation with the initial solution of 
Eq. (|4.7D ; the solution to new equation (see 



Eq. (4.3) ) with the initial condition of Eq. (4.8) 



is plotted in Fig. 11-b, while the solution of 
the BFKL equation with the initial condition 
of Eq. (4.£) is shown in Fig. [l^c. The value of 
Yg is taken to be equal to 3. 

Fig. 0-c 

Figure 11: Solutions to equations as function of k at different energies 



Eq. (4.3) and Eq. (|4.5|) with the initial conditions given by Eq. ([4.7]) ,Eq. 



The solutions to Eq. 

and Eq. ( f4.9| ) are shown in Fig. |ll|. One can see that the solution of the linear BFKL equation (see Fig. 11 



) steeply increases with energy while the solution to the Balitsky-Kovchegov equation shows only a mild 
increase with rapidity ( see Fig. |ll]-a). The solution to the new equation increases with energy but starts 
to fall down at very high energies (see Fig. [lT|-c). From Fig. 11 one can see that the term {N'^)'^ turns out 
to be much smaller than the term AN N'^ in Eq. (p.l|) for small k ( k <2). For large k this term is much 



smaller than the linear term in the equation ( <^ {^'^) )• Therefore, we can conclude a posteriori that 



we can neglect {N'^) -term in Eq. (|2.1|) and reduce this equation to Eq. (4.3) which has been solved. 



The main result of the paper: the estimates for the survival probability ( see Eq. ( 4.24[ )) , is shown in 
Fig. |l^ for two different choice of the dependence of the saturation momentum on b given by Eq. (3.15). The 
first one shown in solid lines, corresponds to Eq. (3.15) where m is chosen the same as in the electromagnetic 
form factor of the proton and value of Qg^ is found from Eq. ( p. 17 ). The second choice was to fix the value 
of to be the same as in Ref. [34] but the value of m is determined by Eq. ( |3.17| ). The survival probability 
is shown by dotted lines in Fig. The qualitative features are seen in Fig. |l^: the value of the survival 
probability is small and its decreases with the growth of energy (rapidity). The smallness stems mostly 
from the steep increase of the solution to the linear BFKL equation while to the considerable decrease 
contribute two factors: the increase of the solution to linear equation and decrease of the solution to the 
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new equation (see Eq. (|4.3| ) ) . The third interesting feature is that the value of the survival probability at 
high energies (large values of rapidity) does not depend on the initial conditions for k >2. It is worthwhile 
mentioning that the errors that stem from the different initial condition turns out to be smaller than 
dS'^/S'^ < 0.1 for any value of k. 




Figure 12: The values of the survival probabilities for leading neutron production in DIS as a function of the photon 
virtuality (q) at different energies (rapidities). The solid and dotted line corresponds to different choice of the value 
for the typical mass in Eq. ( ^.15 ). The solid lines describe the initial condition with m chosen to be the same as in 



the electromagnetic for factor of proton and the value of Qq^ is determined by Eq. ( 3.17 ), while dotted lines present 

is taken to be the same as in Ref. [34] while m is the solution of Eq. (3.17). 



the different choice: Qq^ 



5. Conclusions 



The main result of this paper is Eq. (2T) (see Fig. |^ for the cross section of the inclusive production of 



leading neutrons in DIS. This equation stems from the direct generalization of the approach developed 
in Ref. [18] for the diffractive production in DIS. The asymptotic solution to this equation as well as 
the numerical solution to Eq. ( |2.lD shows that the survival probability defined in Eq. ( |3.lD , is small and 
steeply falls down with energy. We believe that these features of the survival probability are general and 
does not depend on a particular process that we consider. Being the first theoretical attempt to calculate 
the survival probability this paper shows that the survival probability could be as small as 10~^ at high 
energies. 

However, the numbers we need to take with considerable cautions, since these estimates depend cru- 
cially on the assumed impact parameter dependence of both the DIS structure function and the Born 
approximation for the leading neutron inclusive cross section. Unfortunately, the phenomenological analy- 
sis of DIS data based on Balitsky-Kovchegov equation (see Ref. [34] ) was performed neglecting the impact 
parameter dependence. Therefore, to obtain the reliable estimates we need to re-visit the DIS data and 
re-do the analysis using Baitsky-Kovchegov equation the impact parameter depending initial conditions. 
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Therefore, we consider this paper as only the first step to the rehable estimates for the experimentally 
measured cross section. The small value of the survival probability as well as its energy dependence make 
difBcult the task of extraction of the deep inelastic structure function for pion, measuring the spectrum of 
the leading neutron. 
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